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The multifractal theory of turbulence uses a saddle-point evaluation in determining 
the power-law behaviour of structure functions. Without suitable precautions, this could 
lead to the presence of logarithmic corrections, thereby violating known exact relations 
such as the four-fifths law. Using the theory of large deviations applied to the random 
multiplicative model of turbulence and calculating subdominant terms, we explain here 
why such corrections cannot be present. 



1. Introduction 

In fully developed turbulence there is now fairly good evidence for anomalous scaling, 
that is scaling exponents whic h cann ot be predicted by dimensional analysis. Some of this 
evidence is reviewed in lFrischl l|l995|) . This reference also contains a d e tailed presentation 
of the mult ifractal formalism in the formulation of IParisi fc Frischl l)l985|) (henceforth 
PF; see also lBenzi et alJ l)l984j) '). In this formalism, anomalous scaling for structure func- 
tions (moments of velocity increments) is connected by a Legendre transformation to the 
distribution of singularities of the velocity field. An earlier and alternati ve formalism for 
anomalous scaling was introduced by the Russian School of Kolmogorov llObukhovlll962l 
iKolmogorovl 1 1 963: lYaglomlll966|) . In its simplest version it uses a random multiplica- 
tive model for calculating the statistical fluctuations of the ene rgy dissipation on v arious 
scales; the fractal properties of this model wer e disco v ered b y iMandelbrotl l|l974h . The 
bridging of the two formalisms is discussed in iFrischl lll995h in the light of the theory 
of large deviations for the s ums of independ ent identically distributed random variables, 
discovered in the thirties bv lCrarriirl l)l938j) . 

We recall that in its original formulation PF gives an integral representation of the 
structure functions which are then evaluated by the method of steepest descent through 
a saddle point. Specifically, the structure function of order p for a separation I is given 



(1.1) 



We have here used the notation of lFrischl l)l995|) : vq is the r.m.s. velocity fluctuation, £ 
is the integral scale, D{h) is the fractal dimension associated to singularities of scaling 
exponent h and dfi(h) gives the weight of the different exponents For a given p, let us 
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assume that the exponent ph + 3 — D(h) has a minimum ( p , as a function of h and that it 
behaves quadratically nea r this minimum. Standa rd application of Laplace's method of 
steepest descent (see, e.g.. lBender fc Orszadll978j) shows then that, at small separations, 
S p (£) varies as (£/£o)^ p but with a logarithmic prefactor [ln(£/^o)] _1 ^ 2 stemming from the 
Gaussian integration near the minimu m. For p = 3 this log arithmic prefactor is clearly 
inconsistent with the four-fifths law of iKolmogorovl ljl94lj) . one of the very few exact 
results in high-Reynolds number turbulence, which tells us that the third-order (longi- 
tudinal) structure function is given by — (4/5)e£, where e is the mean energy dissipation 

per unit ma ss. 

In iFrischl l|l995|) this difficulty was handled by writing 

tensr^- (L2) 

Indeed, by taking the logarithm of the structure function we change the multiplicative 
logarithmic correction into an additive log-log correction which, after division by ]xi£, 
becomes subdominant as £ — > 0. But if we do not take the logarithm of the structure 
function, is there a logarithmic correction in the leading term whose presence, for p = 3, 
would invalidate the standard multifractal formalism? 

It is also well known that such logarithmic corrections are absent in the random mul- 
tiplicative model, a matter we shall come back to. Actually, as we shall see, the random 
multiplicative model gives the key allowing us to understand why logarithmic corrections 
are unlikely and are definitely ruled out in the third-order structure function. 

In §|21 we recall some basic facts about the random multiplicative model and the way 
it can be reformulated in terms of m ultifract a l sin gularities using large-deviations the- 
ory. This is the theory introduced by ICramerl ( 1938h which allows the estimation of the 
(very small) probability that the sum of a large number of random variables deviates 
strongly from the law of large numbers. f In this way the random multiplicative model 
can be reformulated in standard multifractal language; naive application of PF would 
then suggest the presence of logarithmic corrections. In §|2]we show how a refined version 
of large-deviations theory, which goes beyond leading order, removes the logarithmic cor- 
rections, which are cancelled by other logarithmic corrections in the probability densities. 
In §Q]we summarise our findings, return to the general multifractal formalism beyond 
the specific random multiplicative model and show that the four-fifths law allows us to 
obtain the first subleading correction to the usual multifractal probability. 



2. The random multiplicative model 

In the random multiplicative model (see, e.g.. lFrischlll995l for details) one assumes 
that an integral-scale size cube with side £q is subdivided into 8 first-level cubes of half 
the side, which in turn are divided into 8 2 second- level cubes of side £o2~ 2 , and so on. 
The "local" dissipation at the nth level with scale i = £o2~ n is defined as 

e e =eW 1 W 2 ...W n , (2.1) 

where e is a non-random mean dissipation per unit mass and the Wi's are positive, in- 
dependently and identically distributed random variables of unit mean value. 
The ensemble average (ei) is thus still equal to e, so that the cascade is conservative 

f The reader is only assumed to be fa miliar with a v ery elementary version of large-deviations 
theory, as explained e.g. in Sec. 8.6.4. of IFrischl lfl99Rh : a more detail ed exposi t ion fo r physicists 
intere sted e.g. in the foundations of t hermodynamics can be found in lLanfordl i!973f ) : IVaradhanl 
ljl984f l and lDembo fc Zeitounil ljl998T) are written for the more mathematically minded reader. 
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only in the mean. Since we are interested in describing the incrtial-rangc scaling prop- 
erties (I <C ^o)j we shall mainly focus our attention on high-order generations, i.e. on 
large values of n, where large fluctuations of si are present. As a consequence, the for- 
malism of multiplicative variables leads to the presence of very large fluctuations. The 
correspondence between multifractality and the probabilistic theory is expressed by the 
relationship 

Thus n can be viewed either as the number of W factors determining the local dissipation 
or as the number of cascade steps leading from the "injection" length scale £o to the 
current scale i. 

It is elementary and standard to show that the moments of the local dissipation are 
given by 

, p v - \o g2 (W) 

T p (£) = (e p ) = ((eW 1 ---W n ) p )^s p (W p ) n ^s p (-\ ■ (2.3) 

Then, following the suggestion originally made bvlObukhovl one calculates the 

structure functions at separation £ by the iKolmogorovl l|194l|) expression in which one 
replaces the mean dissipation by its local random value ££, to obtain 

S p (£) « ((ietf/ 3 ) = l p/ % /3 (£) cx & , (2.4) 

where C, p = p/3 — \og 2 (W p ^ 3 ) . Obviously the third-order structure function has exponent 
unity, as required by the four-fifths law and none of the structure functions has any 
multiplicative logarithmic factor. 

There is however an alternative and somewhat roundabout way of evaluating the struc- 
ture functions for the random multiplicative model. First we transform the product of 
positive random variables into a sum by setting Wi = 2~ m \ to obtain 

T v {t) = {{e2- mi ■ ■ • 2- m -)P) = £ p {2- mip ■ ■ ■ 2~ mnP ) 

= e p (2- nxp ) =e p J dxe- nxpla2 P n (x) , (2.5) 

where 

mi + . . . + m„ 



x — 



■_^n (26) 

n 

is the sample mean of n independent and identically distributed random variables and 
P n (x) its probability density function (PDF). When the PDF P(m) of the individual 
variables njj falls off very quickly at large arguments, as is usually assumed in the random 
multiplicativ e model, f all the moments of the miS are finite and the law of large numbers 
l)Fellerll 19681) implies that P n (x) is, f or large n, incr easingly concentrated near the mean 
(m). The theory of large deviations l|Cramerll 19381) tells us roughly that when x ^ (rn) 
its probability falls off exponentially with n as e ns ( x \ where the Cramer function s(x) is 
non-positive up-convex and vanishes at x = (m) . The Cramer function can be expressed 
as the Legendre transform 

s(x) = inf [ax + In Z(a)] (2.7) 

a 

of the characteristic function 

Z{a) ^ J dye- a yp(y) = (e^™) . (2.8) 
f We shall here assume that P(m) falls off faster than exponentially at large \m\. 
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The correct statement of large deviations is that In P n (x)/n tends to s(x) as n — > oo. 
Suppose however we somewhat sloppily write the large-deviations result as 

P n (x) ~ e ns{x) INCORRECT (2.9) 

and use this in 12.5fl . We then obtain an integral representation for T p (£) and thus for 
S p (£) which when evaluated by steepest descent for large n will give not just a power law 
in I but also a multiplicative correction proportional to 1/-%/— ln(£/£o), thus contradicting 

To resolve the paradox we need to extend the large-deviations result beyond the 
leading order. This is called the theory of refined large deviations, first de v eloped by 
iBahadur fc Ranga Raol (^960) and which is reviewed in iDembo fe Zeitounil ljl998j) . In 
the next section we shall show how this can be done by rather elementary application of 
steepest descent. 



3. Refined large- deviations theory and the disappearance of logs 

We now derive the asymptotic expansion for large n of the PDF P n {x) of the sample 
Consider the characteristic function of the sample mean 



Z n {a) = J dxe' aX P n (x) = ( e -"(mi+...+m„)/n) ( 3 J) 
_ f^—ami/n _ _ _ g — am n /n\j _ ^-am/nyi _ ^ — ^ (3 2) 

where Z(a) is the characteristic function for a single variable m, defined in Q2.H)) . Since 
we assumed that P(m) falls off faster than exponentially, Z(a) can be defined for any 
complex a. Hence we can invert the Laplace transform appearing in l|3.1|l by a Fourier 
integral along a contour C running from — ioo to +ioo (see Fig. [2J 

P n {x) = i^- I dfieP'Zof?) . (3.3) 



2i7t J c \ n ; 
We recast (|3.3|l in exponential form 

P n (x) = — [ d /3c (3x+nlnZ( - f3 / n '> = — [ d 7 c"^+ lnZ ^] , (3.4) 
2m J c 2m J c 

with the substitution 7 = (3/n. By (|2.7|l the argument of the exponential has a minimum 
s(x) along the real 7-line at a point a*(x). Taking the contour C through a+(x), the 
argument of the exponential will now h ave a maximum at this point; thus l|3.4|l can be 
evaluated by steepest descent (see, e.g.. lBender fc Orszadll97^) . 
We recall that for an integral of the form 

/(»)= / dyf(y)e n *W ; (3.5) 
Jc 

with a saddle point where <fi' vanishes and where neither / nor <fi" vanish, the large-n 
behaviour is given by 



I +()[ I 



The saddle point formula (|3.ti|l . applied to (|3.4|) gives, after taking a logarithm 



(3.6) 



InPrAx) r ,^lnn ln^TtQ) / 1 
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Figure 1. The integration contour C through the saddle a* in the complex /3-plane. 

where Q > is the second derivative of ax + ln Z(a), evaluated at the saddle point a^x). 
As n does not appear in Q, which is solely a function of x, the right-hand side of i|3.7[l 
is thus structured as an inverse power series in n, except for the first subleading term 
which contains a logarithm. 

Note that expressions such as (I3.7[l are very common in thermodynamic applications 
of large d eviations when dealing with the logarithm of the (very large) number of states 
(see, e.g.. lLanfordlll973|) . 

We can of course rewrite (I3.7|l in exponentiated form as 



In this way we see that P n (x) has a multiplicative \fn correction. Recalling that in the 
random multiplicative model n = — log 2 (€/£o), this correction is just what we need to 
cancel the unwanted logarithms in the structure functions obtained when the incorrect 
form (|2.9(l is used. 

It is important to note that the quantity which goes to a finite limit for large n is 
In P n (x)/n and that for this quantity the correction we have determined is an additive 
subleading term. This is why such terms should be regarded as subleading. Of course the 
cancellation of logarithms for the random multiplicative model cannot take place just 
at the first subleading order, since (|2.4|l is an exact expression and has no logarithms. 
To evaluate refined large deviations to all orders for a general random multiplicative 
model is quite cumbersome and we shall not attempt it here because it would not shed 
additional light on the issue discussed. It can however be done qui te easily for simple 
rando m multiplicative models such as the black-and-white model of iNovikov fc Steward 





(3.8) 




(1964). 
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4. Back to multifractal turbulence 

In multifractal language, the result obtained within the framework of the random 
multiplicative model is that the probability P(£, h) to be within a distance I of the set 
carrying singularities of scaling exponent between h and h + dh is not [I / £o) 3 ~ D ^ h ' dfx(h) 
but is actually given, for small £, by 



which has a subleading logarithmic correction. We recall that it must be qualified "sub- 
leading" because the correct statement of the large-deviations leading-order result in- 
volves the logarithm of the probability divided by the logarithm of the scale. The cor- 
rection is then a subleading additive term. 

It is important to mention that the presence of a square root of a logarithm correction in 
the mu ltifractal probability density has already been proposed bvlMe neveau fc Sreenivasai 



( 1989) on the basis of a normalization requirement; they observed that without such a 
correction the singularity spectrum /(a) comes out wrong; they also po inted out that 
a similar correction has been proposed by Ivan de Water & Schraml l|l988fl in connection 
with the measurement of generalized Renyi dimensions. 

Returning to the multifractal formalism of turbulence, beyond the ran dom multiplica- 
tive m odel, we observe that the usual multifractal ansatz as made in IParisi fc Friscbl 
( 1985) is only about the leading term of the probability, which is easily reinterpreted in 
geometrical language. Hence, it does not allow us to determine logarithmic corrections 
in structure functions. However, if we use Kolmogorov's four- fifths law, we have an addi- 
tional piece of information which implies that the multifractal probability should have a 
subleading logarithmic correction with precisely the form it has in (|4.ip . This improved 
form then rules out subleading logarithmic corrections in any of the structure functions. 

We wish to thank J. Bee, M. Blank, G. Boffetta, A. Celani, A. Dembo, G. Eyink, 
C. Meneveau, K.R. Sreenivasan and A. Vulpiani for illuminating discussions and sugges- 
tions. This research has been partially supported by the MIUR under contract Cofm. 2003 
(prot. 020302002038). 
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